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Abstract. Clustering is a relatively widespread phenomena which takes on many guises across the nuclear 
landscape. Selected topics concerning the study of halo systems and clustering in light, neutron-rich nuclei 
are discussed here through illustrative examples taken from the Be isotopic chain. 

PACS. 27.20. +n 6< A <19 - 21.45. +v Few-body systems 



1 Introduction 

Until relatively recently cluster studies have been confined 
to the region encompassing the line of beta stability where 
the role of a-clustering has long been established Q. As 
[ clustering is expected to manifest itself most strongly near 
. thresholds g, exotic structures might be expected to form 
' in the more weakly bound systems found near the neutron 
, and proton driplines. In practical terms the exploration of 

■ clustering in nuclei far from stability has become techni- 
[ cally feasible over the last 15 years through the advent 
• of radioactive beams and associated techniques. This is 

■ demonstrated most clearly by the discovery and subse- 
quent probing of the nuclear halo. 

In the present paper selected topics concerning the 
study of halo and molecular states in light, neutron-rich 

■ nuclei which have undergone developments since the pre- 
' vious conference in this series are discussed. Following a 

brief review of work on xaiXn molecular states, empha- 
sis is placed on new techniques for probing correlations in 
two-neutron halo nuclei. Very recent work exploring the 
means to produce and detect multineutron clusters is also 
presented. As they exhibit many of the facets of clustering 
and structural evolution far from stability the examples 
discussed here have been selected from the neutron-rich 
Be isotopes. 



2 Nuclear Molecular Clusters 

It has long been established that the a-particle plays 
an important role in the structure of light a-conjugate 
(A=4n) nuclei ||l|. This is a direct consequence of the 
strongly bound character of the '^He nucleus and the weak- 
ness of the a-a interaction, as evidenced by the unbound 
nature of *Be. Whilst an excess of neutrons (or protons) 
may naively be expected to dilute any underlying a-cluster 
structures, theoretical P,W and recent experimental work 



indicate that molecular type structures such as a- 
chains "bound" by valence nucleons may also occur. The 
persistence of such cluster structures in systems lying 
away from the line of beta-stability is well illustrated, as 
will be discussed here, by the beryllium isotopes, for which 
the a-a system may be regarded as the basis. 

From a theoretical point of view, prescriptions such as 
the Molecular-Orbital Model (MO) g or the Two-Centre 
Shell Model 0, in which valence nucleons are added to 
the single-particle orbits arising from the two-centre po- 
tential, provide a successful and conceptually appealing 
framework within which to describe the properties of these 
nuclei. Moreover these orbits may be viewed as the ana- 
logues of the a and 7r-orbitals associated with the co- 
valent binding of atomic molecules. The development of 
fully fledged Antisymmeterised Molecular Dynamics cal- 
culations (AMD), as discussed by Horiuchi in his contri- 
bution to this conference, is of particular interest as the 
A-nucleon system is modelled without any a priori impo- 
sition of an underlying cluster structure. Recent calcula- 
tions, in particular, suggest the existence of two-centred 
structures in the Be, B and C isotopic chains with va- 
lence neutron density distributions exhibiting the features 
of molecular orbitals ||] . 

From an experimental perspective, von Oertzen has 
compiled systematic evidence for the existence of dimers 
in ^~^^Be and ^~^^B. In the case of ^Be, for example, the 
presence of a valence neutron results in a bound (Bor- 
romean) system, the ground and excited states of which 
may be understood in terms of a three-body a:ii:a molecu- 
lar structure. In particular, the rotational bands based on 
the ground and low-lying states exhibit large deformations 
consistent with the associated molecular configurations. 

In the case of ^°Be, the experimental evidence for 
molecular configurations is rather less well documented. 
Beyond the established Oj, 2 J and lj~ - 4j~ states, the 
locations of the J=5 |^ and 6 members of the of the neg- 
ative parity band, as well as the J=4 and 6 members of 
the positive parity band have been postulated following 
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Fig. 1. Spin-energy systematics for states observed in ^ Be 
(from ^). The trajectories for the postulated positive and neg- 
ative parity molecular bands are indicated. 



recent studies of the a-^He breakup of ^°Be* As dis- 

played in Fig. 1 the spin-energy trajectories for the 
bands based on the and 1^ states at ~6 MeV are con- 
sistent with large deformations as expected for molecular- 
like a:2n:a structures. Moreover, the location of the band- 
heads just below the threshold for a + ^He decay is in ac- 
cordance with the considerations of Ikeda describing the 
formation of clusters |^ . Theoretical support for the pos- 
tulated molecular states may be found in AMD calcula- 
tions, whereby well developed Q!:2n:Q; configurations are 
predicted for the and lj~ bands [Q. Coupled-channel 
calculations presented in these proceedings by Ito and col- 
laborators also provide support for the molecular nature 
of the band based on the J . Furthermore, extended MO 
model calculations have demonstrated the that the 
O2 state may be well characterised by valence neutrons in 
occupying the cr-orbital. 

Given the existence of such molecular-type structures 
in ^"Be, the question naturally arises as to the existence 
of similar structures in more neutron-rich systems. In this 
context the dripline nucleus ^^Be has been investigated 
via inelastic scattering at 35 MeV/nucleon. Evidence was 
found in this measurement for the breakup into ^He-|-^He 
and Qf-l-^He of states (J=4, 6, 8) in the excitation energy 
range 10-20 MeV which exhibited spin-energy systematics 
characteristic of a rotational band Moreover the in- 
ferred momenta of inertia (fi,^/23=0.15±0.04 MeV) and 
bandhead energy (10.8±1.8 MeV) of the observed states 
are consistent with the cluster decay of a molecular struc- 
ture which may be associated with a'Aii:a configurations. 
As reported at this conference by Saito et at, an ex- 
periment undertaken very recently at RIKEN using two 
neutron removal from an energetic ^^Be beam has uncov- 
ered evidence for possible 0"^ and 2+ states some 1.7 and 
2.7 MeV above the ^He-l-^He breakup threshold. Interest- 
ingly these states conform well to the spin-energy system- 
atics established in our original study. 



Somewhat indirect evidence for the existence of the 
states observed in the two breakup experiments may be 
found in the observation of highly excited states in the 
^Be(^^N,-'^^N)^^Be reaction which are consistent with the 
spin-energy systematics of a deformed rotational band 
[ p^ . The theoretical investigation of molecular configu- 
rations in ^^Be represents a somewhat more challenging 
venture than the lighter mass Be isotopes. Nevertheless, 
efforts are being made, as evidenced by the contribution 
to this conference by Ito and coworkers and recent papers 
by Itagaki et al. pi| and Descouvemont and Baye [0|. 

As suggested by von Oertzen Q and more recently 
by Itagaki et al Q, the neutron-rich C isotopes may 
be expected to exhibit 3a:Xn cluster structures. In this 
context, we have attempted to observe such states in the 
i2C(i6C,i6C*^i°'i2Be-f reaction at 35 MeV/nucleon 

[ p^ . Careful analysis of the fragment coincidences could, 
however, only put an upper limit of some 30 fjh on the 
yield to states in these decay channels. The inability to 
access such states by inelastic scattering may arise from a 
much smaller overlap between the ^®C ground state and 
the cluster states than in the case of ^^Be. In this con- 
text transfer reactions using projectiles (^'^He) and tar- 
gets (^Li, ^Be, ^^C) exhibiting cluster structure may pro- 
vide a more efficient means of accessing molecular states 
in neutron-rich nuclei. 

Beyond the extension of cluster studies to more exotic 
systems, one of the most pressing issues is the determina- 
tion of partial decay widths. Clearly this is not a trivial 
task, as evidenced by the difficulties encountered in the 
study of clustering in stable nuclei. In the case the light, 
neutron-rich nuclei discussed here, the presence of rela- 
tively few decay channels and the more favourable signal- 
to-background ratio when employing radioactive beams 
may, however, facilitate such measurements. 



3 Correlations in Two-Neutron Halo Nuclei 

Clustering also appears in the guise of neutron haloes in 
ground states near the neutron drip-line. Arguably the 
most intriguing are the Borromean two-neutron halo nu- 
clei (^He, ^^Li and ^''Be), in which the two-body subsys- 
tems are unbound. Such behaviour naturally gives rise to 
the question of the correlations between the constituents. 
Even in the case of the most studied of these nuclei, ^He 
and "'^^Li, little is known in this respect. Over the last three 
years we have developed a number of techniques for ex- 

fring the correlations within two-neutron halo systems 
111- 

We have explored the spatial configuration of the halo 
neutrons at breakup through application of the technique 
of intensity interferometry - an approach first developed 
for stellar interferometry by Hanbury-Brown and Twiss in 
Australia in the 1950's and 60's ||l^ and later extended 
to source size measurements in high energy collisions . 
The principle behind the technique is as follows: when 
identical partiles are emitted in close proximity in space- 
time, the wave function of relative motion is modified by 
the FSI and quantum statistical symmetries — in the 
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case of halo neutrons the overwhelming effect is that of the 
FSI jl^ . Intensity interferometry relates this modification 
to the space-time separation of the particles at emission 
as a function of the four-momenta of the particles through 
the correlation function Cnn, which is defined as, 



(free) 



(n-n) (n-n) + (core-n) 



Cnn(Pl,P2) 



(Pn/dpidp2 
(dn/dpi) {dn/dp2) 



(1) 



where the numerator is the measured two-particle dis- 
tribution and the denominator the product of the indepen- 
dent single-particle distributions As is generally the 
case, the single-particle distribution have been generated 
in our work via event mixing. Importantly, in the case 
of halo neutrons special consideration must be given to 
the strong residual correlations Experimentally care 
needs to be taken to eliminate cross talk . 

As a first step, kinematically complete measurements 
of breakup on a Pb target of ^He, "Li and "Be were 
analysed ||T^. The choice of a high-Z target was made 
to privilege Coulomb induced breakup, whereby the halo 
neutrons may in a first approximation act as spectators 
and for which simultaneous emission may be expected to 
occur. The correlation functions derived from the data, 
assuming simultaneous emission, were compared to an 
analytical formalism based on a Gaussian source pT[ . 
Neutron- neutron separations of r^n^^ — 5.9 ±1.2 fm 
(^He), 6.6 ± 1.5 fm ("Li) and 5.6 ± 1.0 fm (^^Be) were 
thus extracted. These results appear to preclude any 
strong dineutron component in the halo wavefunctions at 
breakup. It is interesting in this context to compare these 
results to the RMS neutron-proton separation of 3.8 fm in 
the deuteron (the only bound two nucleon system). 

The same analysis has been applied to dissociation of 
^■^Be by a C target, in order to investigate the influence of 
the reaction mechanism. A result which hints at a some- 
what larger separation, r™'* = 7.6 ± 1.7 fm, was obtained. 
This raises the question as to whether simultaneous emis- 
sion can be assumed a priori. In principle, the analysis of 
the correlation function in two dimensions, transverse and 
parallel to the total momentum of the pair, would allow 
for the unfolding of the source size and lifetime ||2l[ . Such 
an analysis requires a much larger data set than presently 
available. The two-neutron halo, however, is far less com- 
plex than the systems usually studied via interferometry 
(for example, heavy- ion collisions (l9|). Moreover, the sim- 
ple three-body nature of the system breaking up suggests 
that any delay in the emission of one of the neutrons will 
arise from core-n FSI/resonances in the exit channel, a 
process that may be expected to be enhanced for nuclear 
induced breakup. 

Correlations in three-particle decays are commonly en- 
countered in particle physics and are typically analysed us- 
ing plots of the squared invariant masses of particle pairs 
{M^j,Mf^), with M^y- = {pi+pjY; a technique developed 
by the Australian physicist Richard Dalitz in the early 
1950's In these Dalitz plot representations, FSI or 
resonances lead to a non-uniform population of the sur- 
face within the kinematic boundary defined by energy- 
momentum conservation and the decay energy. In the 
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Fig. 2. Dalitz plot for the simulated decay (see text) of ^^Be. 
In the left panel no FSI are included. 



present case, the core-|-n-|-n system exhibits a distribution 
of decay energies {E^). The associated with each event 
will thus lead to a different kinematic boundary, and the 
resulting plot containing all events cannot be easily inter- 
preted. We have thus introduced a normalised invariant 
mass, 



{mi 



Vfljf 



(rrii + rrij + i?d)^ — {fni + w^)^ 



(2) 



Ei- 



which ranges between and 1 (that is, a relative energy 
, = Mij — mi — mj between and Ed) for all events and 
exhibits a single kinematic boundary. Examples of how n- 
n and core-n FSI may manifest themselves in the Dalitz 
plot for the decay of ^*Be are illustrated in Fig. 2, whereby 
events have been simulated according to the simple inter- 
acting phase-space model described in ref. 1 16 . The inputs 
were an E^ distribution following that measured , the 
Cnn obtained with the C target, and a core-n resonances 
with r = 0.3 MeV at Eq = 0.8 MeV. Note that due to the 
normalisation the (squared) core-neutron invariant mass 
does not present a simple structure directly related to the 
energy of the resonance /FSI . 

As described in ref. [|6| , the data from the dissociation 
of ^''Be on C and Pb targets have been analysed in terms of 
Dalitz plots and compared to the interacting phase-space 
calculations. The results for breakup on the Pb target are 
shown in Fig 3. The good agreement obtained through 
the inclusion of only the n-n FSI suggests that, within the 
presently available statistics, the two-neutrons are emitted 
essentially simultaneously. 

Despite much reduced statistics for the triple coinci- 
dence events needed for this analysis, breakup on the C 
target is consistent with a finite delay between the emis- 
sion of the two neutrons . This may be related to 
the population of core-neutron resonances as suggested 
by more classical analyses of ^^Be-neutron events |2J]. By 
fixing the n-n source size to that derived from the C„„ 
for the Pb target, an average delay in the emission for 
breakup by C of 15011^^; fm/c was estimated. 

Finally, in the context of the influence of the reaction 
mechanism, it is worthwhile noting that whilst some 35% 
of the two-neutron removal cross section on the Pb tar- 
get is attributable to nuclear induced breakup 123], the 
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Fig. 3. Dalitz plot and the projections onto the squared invari- 
ant masses for the dissociation of ^''Be by Pb. The lines are the 
phase-space model simulations with/without (solid/dashed) n- 
n FSI. The inset shows the measured Ed spectrum. 

requirement of two neutrons in coincidence with the ^^Be 
core in the present analysis reduces this to some 15% ~ ap- 
proximately half of the two-neutron removal cross section 
arises from absorption. 



4 Multineutron Clusters 

It is interesting to speculate that multineutron halo nuclei 
and other very neutron-rich systems may contain compo- 
nents of the wavefunction in which the neutrons present a 
relatively compact cluster-like configuration. If this were 
to be the case, then the dissociation of beams of such nu- 
clei may offer a means to produce bound multineutron 
clusters (if they exist) and, more generally study multi- 
neutron correlations. 

To date the majority of searches for multineutron sys- 
tems have relied on very low (typically ~1 nb) cross sec- 
tion double-pion charge exchange (DttCX) and heavy- 
ion transfer reactions (see, for example, refs |^,^). In 
the case of dissociation of an energetic beam of a very 
neutron-rich nucleus, relatively high cross sections (typi- 
cally ~100 mb) are encountered. Thus, even only a small 
component of the wavefunction corresponding to a multi- 
neutron cluster could result in a measurable yield with 
a moderate secondary beam intensity. Furthermore the 
backgrounds encountered in DttCX and heavy-ion trans- 
fer reactions are obviated in direct breakup. 

The difficulty in this approach lies in the direct detec- 
tion of a "^n cluster. One avenue we are exploring is to 
use the detection of the recoiling proton in a liquid scin- 
tillator 1^ . The advantage of a liquid scintallator is that 
neutrons may be discriminated from the 7 and cosmic-ray 
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Fig. 4. Ej,/Er, for A=l (solid line), 2 (dashed), 3 (dotted) 
and 4 (dot-dashed). In the case of A=l, comparison is made 
to single neutron events from ^'^Li -f C. 



backgrounds. Careful source and cosmic-ray calibrations 
permit the charge deposited and hence the energy {Ep) 
of the recoiling proton to be determined. This may be 
compared to the energy derived from the measured time- 
of-flight {En)'- for a single neutron and an ideal detector, 
Ep/En<l; for a realistic detector with finite resolution the 
limit is ~1.4. In the case of a multineutron cluster ("*n) 
Ep can exceed the incident energy per nucleon and Ep/En 
will take on a range of values extending beyond 1.2 — up 
to --3 in the case of A=4 (Fig. 4). 

Analyses of data already available from the breakup of 
^*Be by C show some 7 events with Ep/EnilA (Fig. 5). 
Furthermore all 7 events are in coincidence with ^"^Be frag- 
ments. Intriguingly, preliminary estimates suggest that 
the detection of two or more neutrons in a single mod- 
ule (the most likely source of background) should produce 
much less that one event with Ep/ EnilA in coincidence 
with ^"Be. Further detailed analysis of these events and all 
possible sources of background is presently underway. In 
addition, higher statistics measurements employing beams 
such as ^He are currently being planned. 
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